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Abstract 

Stability of supermassive stars embedded in hot dark matter is 
investigated on the base of the energetic method. Stabihty effect of 
dark mater is compared with rotational stabilization and preference 
of the last factor is advocated. 



1 Introduction 

Supermassive stars, considered hystoricaly as a first model of quasars and 
active galactic nuclei (AGN) suffer from early instability to relativistic col- 
lapse, imposing their short life time (Fowler, 1964; Zeldovich and Novikov, 
1965; Ozernoy, 1966). This instability is connected with a high entropy of the 
matter in such stars, what leads to prevailing of radiation pressure, closeness 
of the adiabatic index to 4/3 and small stability resource. 

A common way to overcome this instability is to consider rotating su- 
perstars, what may postpone the moment of collapse to 3 x 10^ years for 
solid body rotation with angular momentum and mass losses (Bisnovatyi- 
Kogan, Zeldovich and Novikov, 1967), and much longer for a differentially 
rotating star evolving with almost constant angular momentum (Fowler, 
1966; Bisnovatyi-Kogan and Ruzmaikin, 1973). Farther developement of 
a theory and collection of observations led to accreting supermassive black 
hole (Lynden-Bell, 1969) as a most promising model for quasars and AGN. 

*Space Research Institute, Russian Acadamcy of Sciences, Moscow, Russia. Email: 
gkogan@nix .iki.rssi.ru 



1 



Nevertheless, formation of supcrmassive stars on early stages of the Uni- 
verse expansion, their loss of stability with subsequent collapse or explosion 
(Bisnovatyi-Kogan, 1968; Pricke, 1973; Fuller et al, 1986) could be important 
for early formation of heavy elements, observed in spectra of the most distant 
objects with red shift ~ 5, creation of perturbations for large scale structure 
formation, influence on small scale fluctuations of microwave background ra- 
diation (Peebles, 1987; Cen et al, 1993). A necessity of a presence of a dark 
matter in modern cosmological models makes it important to include it into 
stability analysis of supermassive stars. This was done by McLaughlin and 
Fuller (1996), who dealed with nonrotating superstars. Here we consider 
the same problem for rotating superstars, using energetic method, which 
is asymptotically exact for supermassive stars (Fowler, 1964; Zeldovich and 
Novikov, 1967). We show that rotational effects occure to be more important 
for realistic choice of parameters. 

2 Stability analysis 

In supermassivc stars with equation of state P = Pr + Pg = + plZT 
there is Pr 3> Pg due to high entropy of such stars. Besides, such stars are 
fully convective and entropy is uniform over them, so the spatial structure 
is well described by a polytropic distribution, corresponding to 7 = 4/3. In 
the energetic method, considered by Fowler (1964), Zeldovich and Novikov 
(1965), and formalized by Bisnovatyi-Kogan (1966), the spatial profile of the 
star is fixed during perturbations p = pc(l){'m/M), and stability analysis is 
reduced to the set of algebraic equations. When the main terms in total en- 
ergy (without a rest mass) are the energy of radiation and newtonian gravity, 
small terms connecting with plasma energy, general relativistic corrections, 
rotation and a dark matter, are calculated using polytropic, with 7 = 4/3, 
density distribution, corresponding to iscntropic star with pure radiation 
pressure. In this situation influence of a hot dark matter, which density does 
not change during perturbations, should be taken into account by account 
of a newtonian gravitational energy of the star in the dark matter poten- 
tial, because GR effects of a dark matter are of a higher order of magnitude 
(McLaughlin, Fuller, 1996). 

For radiation dominated plasma there is a following expression for the 
adiabatic index, determining the stability to a collapse 
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where P = ^ = In the radiationaly dominated supermassive star there 
is a unique connection between its mass M and entropy per unit mass S 
(Zeldovich and Novikov, 1967) 

where a is a constant of the radiation energy density, and numerical coefficient 
is related to the polytropic density distribution with 7 = 4/3. At the point 
of a loss of stability the critical value of an average adiabatic index < 7 > in 
selfgravitating nonrotating star with account of post-newtonian corrections 
is determined by a relation (Zeldovich and Novikov, 1967) 



3 3 3 3 

Here averaging is done over the stellar volume with a pressure P as a weight 
function, egr and En are post-newtonian and newtonian gravitational ener- 
gies of the superstar, respectively. From comparison between (|ID and (0) we 
get a well known relation for a critical central density of a supermassive star 
stabilized by plasma 

Pc = O.IOtSLm-/^ - 1-8 X f^V''g/cm^ (4) 



Here and below we consider for simplicity a pure hydrogen plasma. New- 
tonian energy of a superstar End in the gravitational field of uniformly dis- 
tributed dark matter with a density pd is written as 

I'M 

Snd = / (fddm. (5) 





The gravitational potential of a uniform dark matter is written as 

- fopy - 1^, (6) 

where Rd is much larger then stellar radius R, and Ma is a total mass of 
the dark matter halo. Stability does not depend on normalization of the 
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gravitational potential so we shall omit the constant value in (|^). It follows 
from d^) and (j^) that during variations End ~ Pc'^^^, while Egr ~ pl^^ and 
~ pI^^ (Zeldovich and Novikov, 1967). Critical value of the adiabatic 
index in the energetic method is obtained in the most convenient way by 
equating to zero a second derivative of the total energy over p^/^. As a result 
for nonrotaing superstar in presence of a dark matter the critical value of an 
average adiabatic index < 7 >crnrot is determined by 

, r ^ , 2 Egr „ ^nd 
< 7 >crnrot=< 7 >crs +ddm = o + q ^ I T' \') 

o o En \^n\ 

Taking into account, that ^at ~ 127r Pr'^dr and (H) - we obtain the 
value of a dark matter correction to the stability criteria 6dm coinsiding with 
the corresponding correction 6dm in the formula (He) from McLaughlin and 
Fuller (1996) (later with subscript mf- (Hgmf))- The relation for a critical 
density in presence of a dark matter is obtained by comparison of (|l]) and 
(i- We get 

= + f = 4.1^1 + -(—) (-) . (8) 

Here relations R = ^^q^26 ' ^^"^ lo^ P''^'^dr = 6.95 x lO'^p^R^, based on 
Emden polytropic distribution with 7 = 4/3, were used (see Bisnovatyi- 
Kogan, 1989). Relation is reduced to {6mf), note a misprinting in the 
sign in {6mf)- 

^■»-^0-'{§-f%r=3.,.10-^{^fp,^,.. (9) 
In Fig.l solution of (H) is presented, as in Fig.2MF- 

3 Stability of rotating stars 

Rotaion was shown to be a main stabilization factor of a superstar. Let us 
consider a rigid rotation, when its energy is a small correction to the energy 
of radiation and the energetic method is a good approach. When losses of an 
angular momentum during evolution are negligible we distinguish between 
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rapidly rotating (RR) and slowly rotating (SR) superstars. In RR case a 
superstar reaches the state of rotational equatorial breaking before loosing its 
dynamical instability, and in SR case instability comes first. If a superstar has 
an angular momentum J, then its rotational energy e^ot ~ 1.25J^p^/^M~^/^, 
and a ratio Srotl^GR remains constant during evolution (Bisnovatyi-Kogan, 
Zeldovich, Novikov, 1967). In presence of rotation and dark matter the 
critical value of the adiabatic index < 7 >crrot is written as 

^ ^ 4 2\eGR\—Srot r,^nd /.^.v 
<l>crrot=- + - , , 2- r, (10) 

and the relations for determination of a critical central density, instead of 
(H), is written as 



As follows from (|11]), a superstar does dot loose its stability when Srot > 
{egrI- This qualitative result, obtained in the post-newtonian approximation, 
remains to be valid in a strong gravitational field and reflects a presence of a 
limiting specific angular momentum aum = GM/ c, so that a black holes with 
a Kerr metric may exist only at a < aum (Misner, Thorne, Wheeler, 1973). 

A RR superstar in a course of the evolution reaches instead a limit of 
a rotational instability, and equatorial mass shedding begins, leading to a 
loss of an angular momentum. Such star will loose the stability when the 
anuglar momentum will become less then the limiting value. The stage of 
a mass loss was examined by (Bisnovatyi-Kogan, Zeldovich, Novikov, 1967), 
where it was shown that this stage may last about 10 times longer, then a 
maximum evolution time to approach the rotational instability point. RR 
star reaches the stage of a rotational instability at different central densities, 
depending on J, but the ratio of rotatonal and Newtonan gravitational energy 
on the mass-shedding curve remains constant (Bisnovatyi-Kogan, Zeldovich, 
Novikov, 1967) Q 

^Note, that in presence of a dark matter halo mass outflow begins at rotational energy 
approximately (1 + bipd/ Pc) times larger, then in (|l^ ) due to additional gravity of a dark 
matter. This correction is small for a considered halo density on the mass-shedding curve 
at J > Jo (see below), and is neglected in this section. 
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erot = 0m725\eN\. (12) 

The energy of a rotating supertar in equilibrium in presence of a hot dark 
matter may be written as 

£eq = —£gas + I^Gr] — ^rot + '^^nd- (13) 

In the main term for a superstar in equihbrium a relation is valid 

£gas = ^\eN\- (14) 

Taking into account (|1^) , (|I^ , we get an expression for an equilibrium energy 
along the mass-shedding curve (with variable J) 

Beg = - (^0.00725 + IsnI + {sgrI + S^w- (15) 

The curve eeq{pc) has a minimum at the central density, determined by a 
relation 

2.8 X 10-3 i^—J p'J' = 3.5 X 10-^ i^-jjj + 5.9 x lO'V + Pd- (16) 

From comparison (|16D and ( p!T]) with account of (|l^) it is clear, that dynamical 
instability cannot occure in the minimum of the mass-shedding curve, and 
after crossing it the evolution proceeds with a substantial mass and angular 
momentum losses. Central density of the superstar in the minimum of the 
mass- shedding curve (|T3p with and without dark matter are represented in 
the fig.l. 

Let us find a parameters of a superstar, at which its critical state is sit- 
uated on the mass-shedding curve. These parameters satisfy sumultanously 
the relations ( |TT1) and (0), what leads to the equation for determination of 
a central density in the form 

2.8 X 10-3 i^—J pT = 3.5 X 10-^ i^-j^j + 12 X 10"^ + Pd- (17) 
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The relation (JT^ is used for determination of an angular momentum of the 
superstar J = Jq with pc from (|16D and J = Ji < Jq with pc from ([T7|) . 
Solution of this equation is also given in fig.l. As may be seen from fig.l, 
the stabilizing effect of rotation on the mass-shedding curve at J = Ji is 
more important, then stabilization by a hot dark matter, which influence 
decreases at increasing of a central density. With account of a longer state of 
the evolution with mass loss until reaching the point of the loss of stability 
at larger J, this conclusion becomes even stronger. 

Determination of the central density at the point of the loss of stability is 
often not enough for characterizing the stabilization. The effect of stabiliza- 
tion is characterized by increase of a life time r of the star until it reaches 
the point of a loss of stability. In turn r depends on the energy loss rate and 
changing of a binding energy of a superstar along its evolutionary path. 



4 Evolution of a superstar 

In massive superstars the temperature and the density at the point of a loss 
of stability are not enough for sufficient nuclear energy production, so grav- 
itational energy is the main source. Luminosity of a superstar is determined 
by its mass and is very close to the Eddington critical luminosity Le = ^^"^^^ 
(Fowler, 1964; Zeldovich and Novikov, 1965) In rotating superstars mass loss 
begins at luminosity less then Le due to the action of a centrifugal force, 
and the presence of a dark matter imply increasing of the luminosity of a 
superstar, because the radiation force must balance the gravity of the dark 
matter in addition to the selfgravity of a superstar. When using Le for esti- 
mation of the evolution time we should take into account possible mass loss 
in the case of a rigid rotation. In the presence of a dark matter we would 
overestimate the lifetme of the superstar by using (??). 

Taking into account the gravity of a dark matter, we obtain that radia- 
tional force balance the total gravity at L = L^e 

i,,.™(l + 54.2^). (18) 

V PcJ 

Nonmonotonial dependence of the binding energy E}, = —Seq on the central 
density may confuse the estimation of a stabilization effect, using only a 
central density data. This is what takes place along the mass-shedding curve. 
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Here maximum of the binding energy corresponds to the central density, 
determined by (16), and at larger central density on this curve in the point 
coinciding with the critical point, determined by (|T7p, the binding energy 
is less, so the evolution time is smaller. Comparison of binding energies in 
these two points is easy do to analytically for small dark matter influence. 
For minimum of mass-shedding curve we get a central density from (O) as 
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At the critical point on this curve it follows from (ITtI) 
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where A = i(0.0075 + f ). The binding energy in the first point, with account 
of dH) is 
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and in the second point, with account of (|20|) is 

/5, 
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We get from comparison of (|21| ) and 



Bbms - e,.s ^ A2(2A-f)3 ^^^^ 



The second term in (|23|) is much smaller then the first, so the minimum point 
with smaller central density (|T8|) has larger binding energy then the point of 
loss of stability (]T9|). Time of the evolution of a superstar until the point of 
the loss of stability r is equal to 

' ^ (24) 

Due to increase of the binding energy with increase of a central density the 
evolution becomes slower with time. In the minimum of the mass- shedding 
curve the central density for M = IO^Mq is about 0.06 g/cm^ (see Fig.l), so 
the increase of the critical luminosity due to the gravity of a dark matter is 
negligible for = 10"^ g/cm^, see (|15). It means, that estimations of the life 
times of the superstar made by Bisnovatyi-Kogan, Zeldovich, Novikov (1967) 
for rotating superstars without a presence of a dark matter remain valid for 
parameters of the dark matter, considered above (McLaughlin, Fuller, 1996). 
The maximal evolution time of a rigidly rotating superstar without a mass 
loss cannot exceed 3 x 10^ years, and the stage of a mass loss until reaching 
the point of a loss of stability could be an order of a magnitude longer. 
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Figure 1: The correction terms Sqr and Sdm + ^GR, the quantities P/6 (line a), 
P/6+{Srot/£N)sh/'i (line b), and 2(3/6+ {Srot/ ^N)sh/'i (line c), as functions of 
the central density of a supermassive star with M = IO^Mq, and dark matter 
density of 10~^ g/cm^. The instabihty points for nonrotating star occure at 
intersection of the correction term curves with the line a. Mass shedding in 
the stable star with angular momentum Jq (see text) occures at intersection 
of correction term curves with the hue b, and critical point on the mass- 
shedding curve is determined by a corresponding intersection with the line 
c. 
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